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1. Introduction 


The Boris algorithm is a popular technique for the numerical time advance 
of charged particles interacting with electric and magnetic fields according to 
the Lorentz force law 0,1 ill- Its popularity stems from simple implemen¬ 
tation, rapid iteration, and excellent long-term numerical fidelity ii- Ex¬ 
cellent long-term behavior strongly suggests the numerical dynamics exhibit 
conservation laws analogous to those governing the continuous Lorentz force 
system Q. Without conserved quantities to constrain the numerical dynamics, 
algorithms typically dissipate or accumulate important observables such as en¬ 
ergy and momentum over long periods of simulated time d- Identification of 
the conservative properties of an algorithm is important for establishing rigor¬ 
ous expectations on the long-term behavior; energy-preserving, symplectic, and 
volume-preserving methods each have particular implications for the qualitative 
numerical behavior SSiSliSIiil. 


A recent Journal of Computational Physics article [12| argued that the Boris 
algorithm is a variational integrator 0 and effectively equivalent to a sym¬ 
plectic integration of Hamilton’s equations for the charged particle dynamics. 
However, this claim contradicts claims made in Refs. which demonstrate 

the volume-preserving property of the Boris scheme. Although symplectic inte¬ 
grators are also volume preserving, volume-preserving integrators are not nec¬ 
essarily symplectic, and Ref. d concludes the Boris algorithm is only volume¬ 
preserving based on direct calculation of a symplecticity condition. Interestingly, 
the calculations in Ref. d do not preclude a discrete variational formulation, 
as claimed in Ref. 12|. The reason for this is that there exists freedom in the 
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identification of a discrete momentum corresponding to a particular state 
with position and velocity v^. The calculations in Ref. Q show that the 
Boris scheme is not symplectic when pfc = mv^ or p^ = mvfc + |A(qfc), where 
m is the particle’s mass, e its charge, c the speed of light and A the magnetic 
vector potential. For variational integrators, the discrete momenta depend on 
the particular choice of discrete Lagrangian, and would likely differ from the 
preceding two definitions. There therefore exists a controversy in the litera¬ 
ture, and identifying the conservation properties of the Boris method requires 
resolving this discrepancy. 

In this letter, we prove that the Boris algorithm does not possess a dis¬ 
crete variational formulation by demonstrating violation of a discrete Helmholtz 
condition 1^. The discrete Helmholtz condition is a necessary and sufficient 


requirement for second-order finite-difference equations to be variational algo¬ 
rithms. The argument for the variational formulation of Boris presented in 
Ref. 12| therefore cannot be correct, and we detail two invalidating flaws in the 
argument. The first flaw involves improper calculation of the discrete Euler- 
Lagrange equations corresponding to the specified discrete action. The second 
flaw is the introduction of truncation error after performing the discrete action 
minimization, as identified as suspect in Ref. |14l |. Ultimately, the Boris algo¬ 
rithm has only been demonstrated to be volume-preserving, and the long-term 
behavior should be attributed to this property. 


2. The Boris algorithm is not a variational integrator 

A popular technique for constructing conservative numerical algorithms is to 
discretize the variational principle underlying the equations of motion [l^ . By 
performing all approximations at the level of the action principle, the numerical 
update inherits several conservation laws of the continuous dynamics, including 
a discrete Noether’s theorem. Centrally, provided certain regularity and con¬ 
sistency conditions are satisfied, performing a time advance using a variational 
integrator is equivalent to using a symplectic integrator on Hamilton’s equations 
[l^ . and is therefore known to exhibit bounded energy error for exponentially 
long times 

Given a finite-difference approximation of an Euler-Lagrange equation, then, 
it is of interest to determine whether the algorithm constitutes a variational in¬ 
tegrator. In the continuous setting, determining whether a second-order differ¬ 
ential equation may be derived as the Euler-Lagrange equations corresponding 
to an action principle is known as “Helmholtz’s inverse problem of the calculus 
of variations” [la, ll7| . Recent work has established a discrete analog to the so¬ 
lution of this problem for determining whether a given finite-difference method 
for second-order differential equations may be derived from a discrete action 
principle, i.e. is a variational integrator [l5|. This section applies the discrete 
Helmholtz condition to the Boris algorithm, and finds that it is only satisfied 
when the magnetic field B(q, f) is a constant. 

A second-order differential equation is specified by a function f(q, v,a, t) 
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and the condition that 




Vt = [0,T] c 


( 1 ) 


Suppose we want to approximate the solution to the differential equation q(t) 
at discrete instances in time tk = 0,/i, ...,Nh, with Nh = T. A second-order 
finite-difference equation is specified by a function P(q, v_, v+, a, t, and the 
condition that: 


P(qfe 


qfc ^qfc-i ^ qfc+i^ qfc ^ - ‘^:^±Bt±^tk,h) = 0, k = 1,N-1, 

( 2 ) 

whereofe is the numerical approximation to the solution q(t = tk)- Following 
Ref. 1^, the finite-difference equation is a variational integrator if and only if 
P satisfies: 


1 /dP dP \ dP dP 

where = (q^, ^ qfc+i-2qfc+qfc-i ^ ^ ^ and the condition must 

hold for arbitrary (qfc,tfc) for all k. Equation ([3]) is known as the discrete 
Helmholtz condition. Note that although the condition in Ref. [l5| presumes a 
one dimensional system for notational clarity, we verified the necessity of this 
condition in higher dimensions. 

The Boris algorithm may be written as [Ij: 


Qfc+i ~ 2qj, -I- qfc_i 


f Qfc-n 


q/c-i 


2h 


X B(qfc) -1- E(qfc). 


( 4 ) 


So, the Boris algorithm is a second-order finite-difference method where the f-th 
component of the vector-valued function P is given by: 

P*(q,v_,v+,a,t,^) = a* - e]i 2 ^ B\q,t) - E\q,t), (5) 

where e is the Levi-Civita symbol, and summation over repeated indices is im¬ 
plied. Evaluation of the terms in the discrete Helmholtz condition Eq. ([3]) yields 
the condition that: 

0 = {B\qk,tk) - B\qk-i,tk-i)) , ( 6 ) 

for all i and j equal to 1, 2 or 3, and for arbitrary (q^-i, q/j, tk-i,tk) in R®. We 
can see, then, that the Boris algorithm admits a variational formulation only in 
the case when B(q, t) is a constant function of position and time. 


3. Critique of the argument in Ref. (l2l | 

Because the Boris algorithm violates the discrete Helmholtz condition, the 
claim in Ref. 12| that the scheme is a variational integrator cannot be correct 
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except in the case of constant magnetic field B(q, t). The first misstep in the 
argument of Ref. [I^ involves mis-calculation of the discrete Euler-Lagrange 
equations corresponding to the specified discrete action. In short, only the first 
component of a two-step map is presented. Although this is a critical flaw for 
the particular argument as presented, we find that the discrete Euler-Lagrange 
equations as presented in Ref. 12| are indeed variational, and we identify a 
different discrete action yielding these discrete Euler-Lagrange equations. The 
second critique is the introduction of truncation error after derivation of the 
discrete Euler-Lagrange equations, and we investigate the claim that the ap¬ 
proximation may be considered equivalent to re-definition of the electric and 
magnetic field quantities. We find that this claim only holds in the case when 
the magnetic field B(q, t) is a constant, consistent with the result of Section [2j 


3.1. Incomplete calculation of the discrete Euler-Lagrange equations 

In constructing a method for non-relativistic charged particle dynamics and 
ultimately arguing for the symplecticity of the Boris update. Ref. 0 approx¬ 
imates a small time interval of the continuous action by specifying a pair of 
discrete Lagrangians. Given a charged particle of mass m and charge e moving 
in electromagnetic fields described by normalized vector potential a = and 
scalar potential ip = Ref. [l^ specifies the following discrete Lagrangians: 


Ld{ql,ql,tk,h) 




^ (a(q^ 4) + a(qL 4 + h/2)) ■ (q^ - q^) - 

’^piql,tk + h/2) (7) 

+ I (a(qfc, + h/2) + a{ql,tk -b h)) ■ (q^ - q^.) - 

^l^{nk,tk + h/2). ( 8 ) 


The corresponding discrete action is given by: 


N-l 

5d(qo,qJ,qo = ■ ■ ■ ,<In) ^ + Ll{ql,ql,tk), (9) 

k=0 

where the dependence on numerical step size h is now implied and the identifica¬ 
tion q^. = q®_|_]^ is made to link between steps. The discrete action is extremized 
by trajectories that satisfy the following discrete Euler-Lagrange equations: 

D2Ldiq!l,ql,tk) + DiLl{ql,ql,tk) = 0 ( 10 ) 

^2i-d(qtqL4)+ E'iEd(q°+i,qUi.4+i) = o. (ii) 


Here, Di is the slot derivative indicating differentiation with respect to the i-th 
argument. These equations correspond to Eqs. (13a, 13b) in Ref. 1^. For the 
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specific discrete Lagrangians defined in Eqs. (171151) . the discrete Euler-Lagrange 
equations become: 


^ - a(q?.,4 + h)) + 

^ Va(q^,4 + /i/2) • (qfc -q?.) - hV^piqj^tk + h/2) = 0 

( 12 ) 

- '^^Jl/2 —2 ~ ^(q^i:^fc+i + ^/2)) + 

i Va(q2,4 + h) ■ (qi^i - q^) = 0. 

(13) 

Equation (23) in Ref. [l^ is presented as the set of discrete Euler-Lagrange 
equations corresponding to the discrete action in Eq. ([5]) . However, c omp arison 
with Eqs. (11211131) reveals that this is incorrect. Equation (23) in Ref. [l2| repre¬ 
sents only the first component of a two-step map, namely Eq. (1T51) or Eq. (ITUl) . 
Trajectories generated by iterating Eq. (1T51) alone do not extremize the discrete 
action given in Eq. dH]), so a discrete variational formulation of the Boris update 
cannot be concluded from the proposed discrete action. 

While the preceding criticism is sufficient to invalidate the specific argu¬ 
ment presented in Ref. [12l] for the discrete variational formulation of the Boris 
scheme, an important question is whether inclusion of the full set of discrete 
Euler-Lagrange equations or selection of a slightly different discrete action could 
repair the argument. Indeed, if one checks the discrete Helmholtz condition on 
Eq. (I12|) . one finds that it is satisfied. A variational formulation for Eq. m 
exists, specifically with a discrete Lagrangian and action given by: 

Ld{qk,qk+i,tk) = 2h ^ (^(qfc) + a(qfc+i)) • (q^+i 

h(p{qk,tk) 

N-1 

Sdiqo, ■■■,qN) = ^ Ld{qk,qk+i,tk)- 

k=0 

To resolve the conflict with the result of Section[2l we turn to the second critique. 

3.2. Truncation of the discrete Euler-Lagrange equations 

Proofs of the structure-preserving properties of variational integrators uti¬ 
lize the discrete variational principle and therefore do not necessarily apply to 
approximations of the discrete Euler-Lagrange equations. Any introduction of 
truncation error after performing the action extremization procedure is liable 
to violate the formal conservation properties, so the variational nature of the 
truncated equations must be independently established. In the case of the Boris 
algorithm, the calculation in Section [2] demonstrates this is not possible. 


- qk) - 

(14) 

(15) 
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In the process of relating the familiar Boris scheme to variational integration 
of charged particle dynamics, Ref. 12| introduces the following approximation: 


1 1 /? (^r\ 

2 -a(qL4 + h)) ss --(q^-q^)-Va(qi,<fc+/i/2)-- —(q^,4+V2). 

( 16 ) 

While Ref. [l^ correctly argues that such a truncation does not affect the local 
order of accuracy of the method, the justification that the approximated equa¬ 
tions retain a variational formulation is insufficient. The Reference identifies 
that justification is required, and states that the approximation may be instead 
interpreted as a re-definition of the electric and magnetic fields. Presumably, 
this means that functions B(q, t), E(q, t) may be defined such that: 


^ B(qfc,4)-hE(qfc,tfe) = ^ (a(qfc_i,4_i) - a(qfc+i,4+i))- 

^Va(qfc,4) • (qfe+i - qfc_i) - V(/3(qfc,4), 

(17) 


for all qfc_i, q/;, q/j+i. However, investigation of the dependence on qfc+i,qfc_i 
reveals that such a definition is possible only in the case of constant magnetic 
field. The left hand side of this equation is linear in the quantities qfc-i,qfe+i 
and independent of so the right hand side must depend on these 

variables in the same manner to be satisfied at all q^-i, q*,, q^+i, tk+i- 

This implies that the magnetic vector potential a(q, t) is linear in its spatial 
dependence and independent in time, and therefore, B(q, t) is constant. Of 
course, the case of constant magnetic field indicates that there are no higher- 
order terms in the Taylor expansion, and that no approximation is performed 
in Eq. (fTBl) . 


4. Conclusion 


In summary, the arguments delivered in Ref. endeavor to establish the 
popular Boris scheme as a variational integrator, contradicting claims to the 
contrary in Refs. @,[ 3 . Because the arguments in Refs [^, [iJI do not rule 
out the existence of a discrete variational formulation for the Boris algorithm, 
we have furnished a proof that no such formulation exists based on a discrete 
Helmholtz condition. The establishment of the discrete variational formulation 
of the Boris scheme in Ref. [13 falls short due to the incomplete specification 
of the discrete Euler-Lagrange equations and (primarily) the introduction of 
truncation error after the action minimization procedure. The same critiques 
pertain to the parallel arguments for the relativistic charged particle algorithm 
in Ref. 12|. The good long-term behavior should then be attributed to the 
volume-preserving properties of the Boris scheme. 
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